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ON SOME RECIPROCAL RELATIONS IN THE THEORY OF NONSTATIONARY FLOWS 

By I. E. Gaeeick 



SUMMARY 

In the theory of nonsfationary flows about airfoils, the 
"indicial lift" function ki(s) of Wagner and the "alter- 
nating lift" function C(k) of Theodorsen have fundamental 
significance. This -paper reports on some interesting 
relations of the nature of Fourier transforms that exist 
between these functions. General problems in transient 
flows about airfoils may be giten a unified broad 
treatment when these functions are employed. Certain 
approximate results also are reported which are of notable 
simplicity, and an analogy with transient electrical flows 
is drawn. 

INTRODUCTION 

There exist at the present time two significant func- 
tions that hare been introduced into the two-dimen- 
sional potential theory of nonuniform motion of airfoils, 
one by Wagner (reference 1) and the other by Theo- 
dorsen (reference 2). "Wagner's function concerns the 
growth of circulation or lift about an airfoil at a small 
fixed angle of attack starting impulsively from rest to a 
uniform velocity v. Theodorsen's function describes the 
lift due to circulation about an airfoil oscillating sinu- 
soidally and moving with uniform velocity v. It is 
the object of this paper to note the usefulness of these 
functions in handling a wide class of problems in tran- 
sient flows about airfoils and to point out certain inter- 
esting relations existing between them. These rela- 
tions are of the nature of Fourier transforms, which 
occur with remarkable abundance in numerous fields of 
physics and which are one of the main studies of the 
recently popular operational-calculus methods. A 
noteworthy analogy between transient hydrodynamic 
Sows and transient electrical flows is also mentioned. 

A third interesting function, which concerns the be- 
havior of an airfoil upon entering a gust, has been 
introduced by Kussner in reference 3, which is an 
excellent survey of the status of the problem of non- 
stationary flows about airfoils. This function and its 
relation to Wagner's function will also be discussed. 

THEORY OF NONSTATIONARY FLOWS 

Wagner's function, &i(s). — Let the chord of the airfoil 
be 2b and let the angle of attack (assumed small) be a. 
Also let the impulsive motion from rest to uniform 
velocity v take place at the origin, *=0 (fig. 1). The 
vertical velocity at the airfoil surface is w=v sine*. 



Then, based on the physical assumption that the ve- 
locity at the trailing edge is finite for all time, Wagner 
derives for the lift, as a function of s=rf/6, 

L=2rbpim:Jc 1 (s) (1) 

The function kj.(s) is illustrated in figure 2. Wagner 
does not derive an explicit analytic expression for 
ki(s) but gives only numerical values. Kussner (ref- 
erence 3) derives a slowly convergent expression in 
series form for jt t (a) that checks Wagner's values. The 
expression is rather long and will not be reproduced 
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FlGDBE 1.— Illustration of nondunenaonal parameter », distance traversed In terms 
of naif chord 6. 

here. It is of great interest to note that the following 
simple expression 

k^s)^!-— (la) 

agrees within 2 percent in the entire range 0<s<o>. 
(Cf. table I.) This expression, which may be regarded 
as a fortunate choice of the author, is especially good 
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Figubb 2.— The function of Wagner. 

in the range s small; it may be reasoned that it repre- 
sents the actual physical state more closely than the 
theoretical solution, since it approaches the steady 
condition (s= °= ) somewhat more slowly. 

It is observed that half the final lift is assumed at 
once and that the lift gradually approaches its asymp- 
tote, 2vbpvw, in agreement with the results for stationary 
flows. 
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The function ki(s) is analogous to what is termed in 
electrical-circuit theory the "indicial-admittance func- 
tion" A{t), which is the current response of a linear 
network to a suddenly applied unit voltage; substitute 
lift for current and unit vertical velocity for unit voltage 
to complete the analogy. 

Theodorsen's function, C(k). — The lift on an airfoil 
oscillating sinusoidally through a small angle of attack 
and moving with uniform velocity v is given by the 
sum of two parts: (a) A classical noncirculatory par£ 
and (b) a part due to circulation. This paper is not 
concerned with the classical part, which consists of the 
virtual inertia terms, the general concepts of which go 
back to the time of Erchoff and Kelvin. 

The steady-state part of the lift due to circulation 
about an oscillating airfoil moving with velocity v has 
been given by Theodorsen as 

L=2irbfr»C(k)Q (2) 

where 

Q=we t <' t is the vertical velocity at the three-quarter 
chord point 1 in complex form, 

and where 

k——> where a is the angular frequency. 
The parameter k permits the substitution of distance s 
in place of time t since ut =< -^^=ks. The function 
C(k) is defined as 



(3) 



TO' 



where Ho <2> and Hi w are Hankel functions (Bessel 
functions of the third kind), or, separated into real and 
imaginary parts and expressed in terms of Bessel 
functions of the first and second kinds, 

C(k)=F(k)+iG(k) (4) 

where 

t Ji(.Ji+Yo) + Y 1 (Y l -J l3 ) 

These functions are illustrated in two ways in figure 3. 

It is important to note that, in the interpretation of 
equation (2), CQc) is considered to operate on the 
function Q. Thus, suppose the actual vertical velocity 
is Wo sinks. This quantity must be axpressed as 
J. P. w a e iU . Then the lift is 

J. P. 2xh|m) <) C(k)e a, 

or 

2irbfmt>o(F t +G*)i sin^s+tan- 1 ^) 

i It is a remarkable fact that the vertical velocity at the three-quarter chord point 
determines the circulation force on the airfoil In oscillator? motions. The lift due to 
circulation acts at the one-quarter chord point. The terms "forward neutral point" 
and "rear neutral point" have been introduced by Kflssner to designate these char- 
acteristto points. 



The lift thus has the same frequency as the vertical 
velocity and both its magnitude and phase are functions 
offr. 

The analogy with alternating currents in electrical 
networks can be mentioned. The function C(k) cor- 
responds to the complex admittance function for alter- 
nating currents (reciprocal of the complex impedance 
function Z(«)). The real and imaginary components 
F and 6 are analogous, respectively, to the alternating- 
current power and quadrature components; the lift duo 
to complex vertical velocity Q corresponds to current 
due to complex voltage E. 
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Figure 3— The function CW-.F-HO ot Theodorsen. 

Superposition principle. — Linearity of the equations 
of Wagner permits one to write, in general, for the lift 
due to a vertical-velocity function tt>(s) (at the three- 
quarter chord point) suddenly applied at the instant 
when 8=0, 

(0)k(s) +£^(8-8^^ (5) 

where w(0) is the value of w at s=0. This result can 
be derived by replacing the function w(js) by a step 
function the envelope of which is 10(3) and by going to 
the limit. The equation is a well-known mathematical 
result and is often employed in electrical-circuit theory 
(reference 4, p. 68). 

Several useful forms of this equation exist, of which 
one is 

-=w(0)^(s)-l- J o /ri(si) ^(s-sOdfi! (5a) 

Equation (5) expresses a noteworthy result, since it 
permits the handling (at least formally or graphically) 
of many problems in transient flows that would other- 
wise be quite laborious. As an example of its applica- 
tion, and as a means of obtaining an interesting result, 
let- ^equation (5a) be applied to the vertical-velocity 
function of the form w(s)=w,fi iU . 
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There results for the lift, when the vertical velocity 
Woe"" is suddenly applied at s=0, 

o^=^i W +vwke a *£k 1 (8 1 )er*«d8 l {5b) 

In order to isolate the transient and steady-state parts, 
the familiar device of separating the interval 0 to 8 
into an interval 0 to °° minus an interval s to «» is used. 

Then, for the steady-state part only (writing s in 
place of Si in the definite integral), 

L=2TbfA-w a ike a ' k^e^'ds 

This equation must agree with Theodorsen's result 
(equation (2)), which may be written in this case 

L=2Tb P w: 0 C(k)e tl ' 

Hence it must follow that 

C{k)=ik^h {8)<r**d8 (6) 

or, in better mathematical form, 

G(k)-l=ikJ^[h(s)-l]e- a 'ds (6a) 

It follows, for the components of C(Jc)=F+iG, that 

=j Q $i(s) — 11 cos&s ds (8) 

These equations can be inverted by the properties of 
Fourier transforms (reference 4, p. 183). Then 



G(k) 
k ' 



sinks dk 



= l+|J a °'^ cosksdk 



0) 
(10) 



Some approximate results. — The integration ex- 
pressed in equation (6) can be performed directly when 

2 

the approximate expression k L (s)^l— is employed. 

Then the following simpler expressions for C{k), F, and 
G, which hold within a few percent, are obtained (cf. 
Jahnke-Emde, "Tables of Functions," p. 80) : 2 



City-ls^ike^Wi-tik) (6b) 
^pS*cos4£ si4&— sin4£ Ci4& 

^jS*cos4£ Ci4&+sin4& si4£ 

The following approximate result is also of interest. 
Let the vertical-velocity function be of the form 
w(s)=w 0 (l— e -7 ') where y determines the rate at which 
w{s) approaches w 0 . Then from (5), using the approxi- 
mate expression for ki{s), 

> The expression (8b) mar be considered the limit, as «-*a>, of the function Cot equa- 
tion (5b)): 

COfc, i) -l=t2it«'*( Ei(-tfi) -E5-(t+») ftH 
The transient term containing t approaches 0 as »-»a>. 



-=l-e-' , -27C- cl+ *^[Ei(4+s)T-Ei(4r)] 



2x&pcuto~ 

Second derivation of equation (6). — The result 
expressed in equation (6) can be demonstrated in still 
another way, which, essentially, is Kussner's treatment. 
(Cf. reference 3, p. 420.) A common artifice in the 
treatment of unit discontinuities is to represent the 
unit "jump" function 

l(s)=lfors>0 

l(s)=0for s<0 
by the following integral in the complex k plane, which 
can be evaluated by residue theory: 

Here the hook integral means integrate from — » to 
+ 00 . bypassing the singular point at the origin by a small 
semicircle from below. In effect, a spectrum analysis 
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Piotrax 4.— The function ii(«) of Knssner. 

of the unit jump function has been made as a limit of 
a sum of exponential terms. To each exponential term 
of vertical velocity there corresponds the lift given by 
multiplication with Theodorsen's function 0(k). Then, 3 
by addition 

"2«'-Cr+ 

G{k)-l 



dk 



=0 [s<0] 



dk [s>0] 



Separating C{k) into F{k)+iG{k) and noting that 
G{—k)=F{k) —iG(k), there results 



_2 r 



0 k 

F(k) 



cosks dk 



sin&s dk 



which check equations (9) and (10) and hence also 
equations (7) and (8). 

Kussner's function, £*(s). — Kussner has derived a 
function k 2 {s) (fig. 4), which gives the lift on an airfoil 
as it 'penetrates into a sudden vertical-gust region with- 
out change in direction. If the change in vertical 



> The relation between Cfjt) and ti(t) Is expressed in the 
CfjfcHW — ki(t), where ik is the operator ifit and the e 
function C(k) operating on the unit Jump function 1(f).' 



itional calculus as 
Is Interpreted "the 
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velocity in the gust region is v> (assumed constant), the 
change in liftrias the airfoil leading edge penetrates the 
gust region is 

L=2-rbpvwk s (8) (11) 

Note that the function & a (s) is denned for fi> — 1, and is 
0 for s<— 1, i. e., before the leading edge has pene- 
trated the gust region. 

In order to obtain the relation of k 2 (s) to &i(s), note 
that the operational equivalent of & 2 (a) (in the same 
sense* that C(k)l(8)=h(s)) is CQc) [J a Qc)+iJi<Vi\ 
(reference 3, p. 420). This latter function describes 
the steady-state lift due to circulation on an airfoil 
moving with uniform velocity v and whose vertical 
velocity is oscillating sinusoidally but progressing in 
the form of waves from point to point, i. e., the vertical 
velocity is of the form 

where x defines any point of the airfoil measured from 
the center. 

Then, the following relations hold, writing a for 8+1, 
CQc) l(«r) 

[</«,(*) +*«/,(*)] lM=S(<7)=|+iarcsin (,-1) 



+7VI- ('-!>*" [6<<<2] 
= 1 [<r>2] 

The second relation is given by Nielsen in "Handbuch 
der Oylinderfmictionen," page 197. By superposition, 
there results for the combined operator 

CQc) M><*)+tJ,(*)] l(o-), 



where 



Hence, 



k 3 (s) =£k 1 (a- \)S' (\)d\ 
=0 [X>2] 



A " a(s)= rX' Ai(ff ~ x) V^ dx [ ° <<r<2] (i2) 



/2-X 



d\ [o->2] 



or, expressed in terms of s, 

^(«)=^/_'*l(8-Sl) A /^(fol hl«l] 

== iX! 1 * i(s - Si) Virl dSi [s>1] (12a) 

The effect of an arbitrary gust function w{<r) can be 
written directly by superposition (o-=s-j-l) 

2^=w(0)ft«(o-) +J\(<r l )^w( ff -<rJd<T 1 (13) 



« For a correction to an error in sign that exists in reference 3, p. 420, consult refer- 
ence 5. The values of k»(s) here given therefore need to be modified; the corrected 
values, inoluding the apparent-mass effect due to change of shape, are presented in 
reference S. 



2 

The approximate expression ki(s)£il—£^ may bo put 
into equation (12). Then (cf. table I), 



*S,W«2^±f-I l»2] 

e^V^-cO-cos-'a-o-) 
+4- 



(12b) 




CONCLUDING REMARKS 

It has been shown that- the functions CQc) and &i(s) 
are of considerable significance in the theory of non- 
stationary flows. To a certain extent, the results arc 
formal since many of the analytic properties of these 
functions he hidden in their complicated structures. 
Further mathematical studies of the function CQc) 
when k is a complex variable would appear to be desir- 
able, since this function is associated with the lift due 
to a general damped sinusoidal motion of the airfoil. In 
many problems in transient flows, it is therefore of 
value to employ the approximate expressions for 
CQc), ki(js), and &»(«) given and thus to obtain quickly 
a simpler perspective of the problem. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., March 28, 19S8. 
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TABLE I.— VALUES OF *,(«) AND *,<») 





'*»(»> 




»*>(•> 


•M«)«j>ix.t. 


£-1 


0 


0 


0 


0 


—.5 


0 


0 


.328 


.328 


0 


.5000 


.5000 


.408 


.468 


.6 


.6567 


.5556 


.508 


.508 


1 


.6006 


.0000 


.031 


.033 


2 


.6693 


.6067 


.003 


.090 


3 


.7195 


.7143 


.737 


.732 


4 


.7582 


.7500 


.772 


.764 


6 


.7880 


.7778 


.799 


.789 


6 


.8125 


.8000 


.822 


.809 


7 


.8325 


.8182 


.810 


.820 


8 


.8486 


.8333 


.855 


.810 


. 8 


.8626 


.8161 


.80S 


.862 


10 


.8715 


.8572 


.877 


.800 


20 


.9321 


.9107 


.934 


.918 


CO 


L0000 


LOO0O 


1.000 


1.000 



i Values of ii(i> taken from references 1 and 3. 
1 From equation (la). 

> Values of *,(») calculated from equation (12). 
4 From equation (12J>). 



